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Abstract: A general criterion is formulated in order to decide if recycling of B-cells exists 
in GC reactions. The criterion is independent of the selection and affinity maturation 
process and based solely on total centroblast population arguments. An experimental test 
is proposed to verify whether the criterion is fulfilled. 
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1 Introduction 



The affinity maturation process in germinal center (GC) reactions has been well charac- 
terized in the last decade. Despite a lot of progress concerning the morphology of GCs 
and the stages of the selection process, a fundamental question remains unsolved: Does 
recycling exist? Recycling means a back differentiation of antibody presenting centrocytes 
(which undergo the selection process and interact with antigen fragments or T-cells) to 
centroblasts (which do not present antibodies and therefore do not interact with antigen 
fragments but proliferate and mutate). 

The existence of recycling in the GC has important consequences for the structure of 
the affinity maturation process in GC reactions. The centroblasts proliferate and mutate 
with high rates in the environment of follicular dendritic cells. During the GC reaction 
they differentiate to antibody presenting centrocytes which may then be selected by inter- 
action with antigen fragments and T-cells. If the positively selected centrocytes recycle to 
(proliferating and mutating) centroblasts the antibody-optimization process in GCs may 
be compatible with random mutations of the centroblasts. On the contrary, one-pass GCs 
(without recycling) seem to be inconsistent with an optimization process which is based 
on undirected mutations. As a consequence, the search for mechanisms which facilitate 
the finding of the optimal antibody-type in GC reactions through somatic hypermutations 
is directly correlated to the existence of a centrocyte recycling process. 

The recycling hypothesis was ffist formulated in 1993 [Q] and there exist several ap- 
proaches to verify it. (e.g. |^, Q). Experimentally, recycling of B-cells has never been 
observed directly. This may be due to a small number of recycled cells compared to the 
huge total number of cells in the GC. An indirect experimental evidence for cell recy- 
cling was found by observing the rates of centrocyte apoptosis and the total cell population 
decrease after injection of a second antigen (different from the antigen which originally ini- 
tiated the GC reaction) into a GC which already was in a later phase of its reaction. Both, 
apoptosis of centrocytes as well as the total cell population decrease were faster than in 
unperturbed GC reactions. Interestingly, the total cell population decrease was too fast to 
be solely explained by centrocyte apoptosis. This observation lead to the conclusion that 
a feedback mechanism from centrocytes to centroblasts should exist. Recycling is a good 
candidate for such a mechanism. 

The major conceptual problem in the theoretical contributions is the number of as- 
sumptions which are necessary to describe the affinity maturation process. The recycling 
hypothesis was formulated because it seemed unlikely that a complex optimization process 
(involving up to about 9 somatic hyper-mutations p, |]) may occur in a one-pass GC, 
i.e. without recycling: Starting from about 12000 cells which mutate in a high- dimensional 
(shape) space the probability (mutation is likely to occur stochastically 0, ||) for at least 
one cell to find the right mutant is very small. It would help, if the already affinity matu- 
rated cells could reenter the proliferation process in order to multiply and to enhance the 
chance to find the optimal mutant. 

Due to the origin of the recycling hypothesis, the efforts to show its existence were based 
on the influence of recycling on the maturation process. Consequently, it was found that the 
number of maturated cells becomes substantially larger with recycling than without |^ . A 
quantitative description showed, that in order to bring the affinity maturation process into 
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accordance with experimental observations, the recychng probabihty of positively selected 
centrocytes should have as high values as 80% Q. In both models important assumptions 
were made to represent the antibody-types in a shape space (affinity classes, dimension of 
shape space) and to allow a successive affinity maturation (affinity neighborhood). 

Here, a new and simpler perspective is presented. Not the influence of recycling on 
the maturation process is examined but its influence on the total centroblast population. 
The centroblasts are believed not to present antibodies, so that they do not participate at 
the selection process via interaction with the antigens presented by the follicular dendritic 
cells. Therefore, in GCs without recycling the total centroblast population is independent 
of affinity maturation and selection. However, this statement may depend from the sys- 
tem under consideration. For instance, the centroblasts in GCs during splenic immune 
responses of mice has been observed to present antibodies 0. 

The analysis starts from the presupposition that recycling does not exist. With a 
minimal set of assumptions (see Sec. ^ the implications for the time course of the total 
centroblast population are described in Sec. |^. The dependency of the results on some 
characteristic properties of various GCs is analyzed (see Sec. |3.1| ) and an experimental test 
of the implications is proposed (see Sec. ||). 

2 Requirements 

2.1 General assumptions 

In this Gedanken experiment the total centroblast population B{t) in a one-pass GC 
(without recycling) is considered. Total means, that the encoded antibody-type is ignored 
and all centroblasts enter the total population with equal weight (despite any differences 
in affinity to a presented antigen). 

Presuming a one-pass GC no centrocytes are recycled to centroblasts, so that the total 
centroblast population becomes independent of the number of centrocytes. Basically, only 
centrocytes present antibodies and thus may undergo a selection process. As recoil effects 
to the total centroblast population are excluded, this magnitude is also independent of 
any properties or dynamics of the selection process. Especially, no influence of the already 
achieved affinity maturation during the optimization process exists. Also, neither the 
number of positively selected centrocytes, nor the speed of the selection process, nor the 
interaction with T-cells, nor the fate of the rejected centrocytes, nor the production of 
plasma and memory cells have any influence on the total centroblast population. 

The total centroblast population changes exclusively according to centroblast prolif- 
eration and the differentiation of centroblasts to centrocytes. Therefore, B{t) follows a 
differential equation of the following type 

^(t) = -fimt) = p{t)B{t) - g{t)B{t) , (1) 

with time t, centroblast proliferation rate p{t), rate of centroblast differentiation to cen- 
trocytes g{t), and total rate function f{t) = g{t) —p{t). 

Cell proliferation has unquestionably to be described by an exponential increase of the 
population. Therefore, a linear proliferation term is used at every time of the GC reaction. 
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In the same way the differentiation term should be directly proportional to the centroblast 
population at every moment. Therefore, nonlinear terms do not enter the total centroblast 
population dynamics Eq. (|1|). 

If the centroblasts participate at the selection process as observed in splenic immune 
responses in mice 0, i.e. if the centroblasts may, as the centrocytes do, undergoe apoptosis, 
Eq. Ill remains correct but with a slightly modified interpretation: The growth of the number 
of centroblasts is inhibited by their apoptosis. Therefore, the function p{t) denotes the 
prohferation rate of centroblasts reduced by the centroblast apoptosis rate. This should 
be beared in mind for the interpretation of the results. 

Note, that no place is determined where proliferation or differentiation occurs. There- 
fore, the following argument is not restricted by the appearance or not- appearance of dark 



and light zones |T0|, M during a GC reaction. 



2.2 About the rates 

It is assumed that after an initial pure proliferation phase of some seeder cells, the total 
centroblast population does not further increase. This is in accordance with measurements 
of the follicle center volume, which show a peak 3 or 4 days after immunization [|TT|. As a 
consequence, the function f(t) in Eq. (|lD is positive during the whole GC reaction (started 
at t = 0), i.e. after the first pure proliferation phase. 

The proliferation rate is probably constant during the whole GC reaction. This has 
not been shown experimentally in all details but is well established by a lot of independent 
data taken at different moments of the GC reaction [^, jl^, 0. All these data lead to the 
same large proliferation rate of 

^.J-. (2) 
ln(2) 6h ^ ^ 

Nevertheless, for the following considerations the value of the proliferation rate will not be 
fixed. In order to assure generality of the argument proliferation rates according to Eq. 
will be discussed as special case only. Therefore, the results of the following analysis apply 
also to GCs with a centroblast proliferation rate that varies in time. 

The situation is more complicated in the case of the differentiation rate g{t). This value 
is not known experimentally, such that its absolute number as well as its time course are 
unclear. From the injection of BrdUrd (5-bromo-2'-deoxyuridine) it is possible to identify 
the cells which were in cell cycle during the preceding 6 hours. In this way it was found 
in GCs with established dark and light zones, that the centrocyte population is renewed 
from centroblasts every 7 hours |TT|. It follows for the differentiation rate that 



ln(2) ^ 7h- 

Because of the lack of more detailed information, different possibilities for the time de- 
pendence of the centroblasts to centrocytes differentiation rate will be considered in the 
following argument. 
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2.3 A standard germinal center 



Before analyzing the implications of this Gedanken experiment a typical GC has to be 
defined as reference system. The GC reaction is initiated by a few seeder cells in the 
environment of follicular dendritic cells (FDC) |1^, |I^. These cells proliferate for 
about 3 days, and give rise to about 12000 centroblasts 0]. At this time to = Oh the 
intrinsic GC reaction starts, including a continuation of proliferation and, additionally, the 
differentiation of centroblasts to centrocytes. The analysis made here does not consider 
the first phase of pure proliferation, but deals solely with the GC properties in the second 
phase. The life time T of a germinal center is about 18 days (T = 432 h) ||Tl], |1^, in which 



the first phase of pure proliferation is not included. At the end of the GC reaction, only a 
few cells remain in the FDC environment [|TT], 0, which is assumed to be about 5 cells in 
a standard GC. The sensitivity of the results on these assumptions will be analyzed (see 
Sec. O 



3 Analysis of the centroblast population 

The different time courses considered for the total centroblast population of a standard GC 
(defined in Sec. p.3|) are shown in Fig.|l|. An exponential decrease corresponds to a constant 
rate function (see Fig.|^ dashed line). A linear population decrease is achieved with the 
rate function Eq. (|T5D (see Fig. |] dotted line) which is derived in the Appendix. Finally, 
using a quadratic rate function Eq. ( p!7D (see Fig. |^ full line) the centroblast population 
remains unchanged for the major duration of the GC reaction and is steeply reduced at 
the very end. 

Experimentally, the time course of the GC total centroblast population seems to lie in 
between the exponential and the linear scenario [0, |l^. Anyhow, the scenario with a long- 



lasting constant centroblast population, i.e. with the quadratic rate function, is certainly 
a reasonable upper bound for possible time courses of the total centroblast population 
(Fig.[^ full line). On the other hand, a population decrease that is faster than exponential 
is unlikely. 

It follows that during the major life time of GCs (t < 15 d) the condition 

fit) = git) -Pit) < 0.018//1 (4) 

holds (see Fig.^. On the other hand the requirement of a monotonically decreasing 
centroblast population leads to the bound fit) > 0, i.e. git) > pit). Taking these results 
together it follows 

Pit) < git) < p(t) + 0.018//i, (5) 

which is a very powerful condition. This means that in an experiment the centroblast to 
centrocyte differentiation rate git) has to fulfill this condition during the first 15 days of 
the GC reaction. If it turns out that the above condition is not guaranteed experimentally, 
then - assuming the rather weak requirements made to be valid - a GC reaction without 
recycling of centrocytes to centroblasts is ruled out. 

Note, that assuming the proliferation rate to adopt the value Eq. (|^), the above condi- 
tion can be reformulated. It follows a lower bound for the differentiation rate which is in 
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f=const/h 

f=-V(t-c), 

f=lh/(t-c)^ 
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Figure 1: Considered time courses of the centroblast population during GC reactions. The 
model results remain valid for all GCs with centroblast populations within these extreme sce- 
narios. 



accordance with Eq. (|^): 

1 git) 1 

6h ln(2) 5.19/1 ^ ^ 

This means that if GCs without recychng exist, a centroblast shall (in average) take not 
less than 5.19 hours to differentiate into centrocytes. 

Supposing that the linear centroblast population decrease is the realistic scenario, the 
upper bound becomes even more powerful. During the first 10 days the condition 

p{t) < g{t) < p{t) + 0.005//1 (7) 

must hold in a GC reaction without recycling. With the proliferation rate in Eq. this 
corresponds to an average centroblast differentiation time of at least 5.75 hours. This result 
demonstrates, that the condition for a one-pass GC may become stronger than in Eqs. 
and (@), depending on the characteristics of the total centroblast population decrease. 
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f=const/h 

f=-V(t-c,), 

f=lh/(t-C2)' 
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Figure 2: The function f{t) = g{t) — p{t) corresponding to the different time courses of the 
centroblast population in Fig.|I]. The GC life time is T = 432 /i. During the major life time of 
the GC the rate function respects the upper bound f(t) < 0.018//i (see Eq. 



3.1 Variation of the GC properties 
3.1.1 The case of constant rates 

In order to check the robustness of the results Eqs. (H) and (|) the most critical scenario 
(i.e. with weakest conditions) is considered in some more detail: the scenario with constant 
rates. If the proliferation and differentiation rate are constant, the centroblast population 
model reduces to a pure exponential decrease of the population 

B{t) = B{0) expiip-g)t) , (8) 

where the constant rate g is unknown. Focusing on the population after T = 18 days of 
the GC reaction, the rate function 

is found. The resulting values for the differentiation rate are shown in Fig. ^ in dependence 
of the GC life time and for different numbers of remaining cells at the end of the reaction. 



7 



As expected the rate becomes larger for shorter GC hfe times T. Even supposing a rather 
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Figure 3: The values of the (constant) rate function / in the case of an exponentially decreasing 
centroblast population. The dependence of condition Eq. ^ of the GC life time and of the 
number of remaining centroblasts at the end of the GC reaction is shown to be weak (see also 
Eq. |10| for the most pessimistic case). 

short hfe time of 12 days and further supposing that only one cell remains at the end of 
the GC reaction, one gets an upper bound of 



p{t) < g{t) < p(t) + 0.033//i 



(10) 



for possible values of the differentiation rate. Incorporating the proliferation rate Eq. 
it follows 

J_ ^ 9{t) 1 



< 



11) 



ln(2) ~ 4.7/1 

As a consequence, the statement that in one-pass GCs (without recycling) the centroblast 
to centrocyte differentiation rate should respect the above bounds Eqs. (|^) and (^ is not 
altered dramatically by variation of the GC properties. 

The initial number of centroblasts 5(0) was not varied so far. Note that only the 
ratio B{T)/B{0) enters Eq. (P), so that a variation of the initial number of centroblast is 
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equivalent to the variation of the final number of remaining cells B{T). Therefore, the 
constraints Eqs. ([T0|) and (0) are also valid with e.g. 60000 initial and 5 remaining cells. 
Lower numbers of initial cells allow even stronger bounds than the above ones. 

3.1.2 The case of linear centroblast decrease 
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Figure 4: The time course of the rate function / = g—^ (corresponding to a linear centroblast 
population decrease) for different GC life times T and for different numbers of remaining 
centroblasts -B(T) at the end of the GC reaction. Condition Eq. ^ turns out not to depend on 
the number of remaining centroblasts. The rate function is sensitive to the total GC life time 
in the late phase of the reaction. However, condition Eq. | remains valid for practically the 
whole GC life time. 



In the case of a linear centroblast population decrease (see Fig.|T] pointed line) the 
corresponding rate function jif) (see Eq. (^)) is analyzed for various GC life times T 
and numbers of remaining blasts at the end of the reaction B{T). The result is shown in 
Fig. ^. Condition Eq. (|^) is not affected at all by a variation of the number of remaining 
cells at the end of the GC reaction in a range of three orders of magnitude (the values 
chosen for B{T) are 0.1, 5, and 100). The three corresponding curves are practically not 
distinguishable (see Fig. |). 
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The situation is different for tlie variation of tlie GC life time T. Here, the period in 
which condition Eq. (^) remains vahd is prolonged for long living GCs and shortened for 
short living GCs (see Fig.|^ dashed and pointed line, respectively). For a very short life 
time of 12 days (plus 3 days of pure proliferation) the condition remains valid for 4 days 
only (7 days after immunization). However, condition Eq. (|^) remains valid even in this 
case for 11 days, so for approximately the whole GC life time. 



4 Conclusions 

The above analysis is summarized with the statement that in a one-pass GC (without recy- 
cling) the centroblast to centrocyte differentiation rate should respect the strong condition 
Eq. (!]) p(t) < g(t) < pit) + 0.018//i during the first 15 days of the GC reaction (with 
standard properties which are defined in|2]^). For a standard proliferation rate Eq. (Q) of 
ln(2)/6/i this translates into an average time for centroblasts to differentiate into centro- 
cytes which cannot be substantially shorter than 5 hours. If this condition is not fulfilled, 
recycling has necessarily to be present in GC reactions. Note that, the other way round, if 
the average differentiation of centroblasts to centrocytes takes more than 5 hours, recycling 
is not necessarily absent. 

This statement is based essentially on two (weak) assumptions: Firstly, linear pop- 
ulation dynamics (see Eq. (0)) for the total centroblast population in GCs. Secondly, a 
monotonic decrease of the total centroblast population. The conclusions apply universally 
to GCs with various characteristics due to the independence of all speculations on shape 
spaces, selection processes, and affinity maturation. However, one should be aware that 
still some assumptions may be questioned. Remind, for instance, that for GCs with se- 
lectable centroblasts (that may undergoe apoptosis) the rate p has to be interpreted as the 
rate of centroblast proliferation reduced by the centroblast apoptosis rate. 

In order to decide if recycling does exist in GC reactions one has to determine the 
centroblast to centrocyte differentiation rate experimentally at least at one representative 
moment of the GC reaction. The moment has to be chosen such that on one hand the 
GC reaction is already well established. On the other hand, the later the measurement is 
done the less is known on the proliferation rate. In addition, one approaches the regime 
where the above condition looses its validity (see Fig. Q). Therefore, an optimal time will 
be about 8 days after immunization. However, if one expects the centroblast proliferation 
rate to vary substantially during this central phase of GC reaction, the standard values 
for p cannot be used. In this case, one would have to determine not only the centroblast 
differentiation rate but also the centroblast proliferation rate. 

Technically, the centroblast differentiation rate may be measured by labeling studies of 



B-cells in the GC which are in cell cycle as it was done in |TT . 

Alternatively, a systematic and time resolved analysis of the total centroblast population at 
different moments of the GC reaction may improve the criterion (see e.g. Eq. (|^)) because 
one could distinguish which of the cell population scenarios presented in Fig. |l] corresponds 
to real GC reactions. This may be done by the measurement of the GC volume at several 



times of the GC reaction as in |11, 17, ^ 
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Appendix 



The differential equation Eq. (|l|) is formally solved by 

B{t) = B{to = Oh) exp(^- dxf{x)^ . (12) 

In the following, the border conditions (life time of GC T, and centroblast population at 
to = Oh and T) are incorporated into the rate function /. 

To achieve a linear decrease of the centroblast population, the second time derivative 
of the population is required to vanish. From Eq. (|l|) and using Eq. (|12D it follows 

(fB_ _ _d£^ _ dB_ 
'd^ ~ di ~dt 

= -f 5(0) exp J^dxfix)^ + fB{0)f exp J^dxf{x)^ . (13) 
This expression vanishes for 

f , (14) 



dt 

which is solved by 

fit) = -r^ • (15) 

t — Ci 

The integration constant ci is fixed by the condition that the final population at time T 
is given by B{T). With Eq. (|T2|) this leads to 

T 

^1 = . B(T) ' (16) 
3(0] 

which is a positive constant for B{T) < B{0). 
In the case of a quadratic rate function 

[t - C2Y 

an analogous calculation which incorporates the final centroblast populations leads to 



T 

C2 = + 



2 \ 4 Inf:^ 



Other solutions are formally possible but do not lead to biologically reasonable centroblast 
populations. 
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